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(f) ' Abstract. The two-pion production reaction pp — -> ppix + -K was measured with a polarized proton beam at 

T p w 750 and 800 MeV using the short version of the COSY-TOF spectrometer. The implementation of a 
delayed pulse technique for Quirl and central calorimeter provided positive ir + identification in addition to 
the standard particle identification, energy determination as well as time-of-flight and angle measurements. 
| Thus all four-momenta of the emerging particles could be determined with 1-4 overconstraints. Total 

and differential cross sections as well as angular distributions of the vector analyzing power have been 
obtained. They are compared to previous data and theoretical calculations. In contrast to predictions we 
find significant analyzing power values up to A y = 0.3. 
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1 Introduction 

The Z\(1232) resonance is the most dominant resonance 
in the pion-nucleon system. Therefore pion production in 
nucleon-nucleon collisions is strongly affected by its ex- 
citation if kinematically and dynamically allowed [HE]- 
This is not the case in two-pion production near threshold. 
Since the A resonance decays by emission of only one pion, 
single- A excitation can contribute to two-pion production 
only if a second pion is produced associatedly. However, 
near the two-pion production threshold it is very unlikely 
that the second nucleon involved in the collision process 
gets excited, too. Hence the associated production of the 
second pion can only proceed by non-resonant rescatter- 
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ing, which however, is a very weak process as calculated 
by Ref. 0. 

This situation provides a unique opportunity to study 
excitation and decay of the Roper resonance AT* (1440), 
because this process provides the only possibility for reso- 
nant two-pion production at energies close to the two-pion 
threshold - as has been predicted in detailed theoretical 
calculations p] and as has been demonstrated by the first 
exclusive measurements performed by the PROMICE/ 
WASA collaboration at the CELSIUS ring in the energy 
range of T p = 650 - 775 MeV [4ll5ll6j . 

The measured angular distributions in the center-of- 
mass system (cms) are all consistent with isotropy with 
the exception of the proton angular distribution, which is 
shaped by the dominating a exchange between the two 
colliding protons - as demonstrated in Ref. [4] and pre- 
dicted by Ref. [3]. From this we learn that except for the 
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Fig. 1. Sketch of the short COSY-TOF setup used in this experiment showing (a) the full detector arrangement with inserts 
for (b) the start detector region with start wedges and hodoscope displaying both side (left) and front (right) views, (c) the 
central hodoscope ("Quirl") as the stop detector for TOF measurements and (d) the central calorimeter. The positions of the 
subdetector systems (b) - (d) in the full setup are indicated in (a). For a detailed description see text. 
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relative motion of the two protons all ejectiles are in 
relative s-waves to each other at beam energies close to 
threshold. This is not surprising, since the leading mech- 
anism for two-pion production, which is the excitation of 
the Roper resonance and its subsequent decay iV* — > Na 
into the a channel (itk)i=i=o, produces only particles in 
relative s-waves. Also the Roper decay via the A resonance 
N* — > Ak — > Nir + ir~ produces predominantly only rel- 
ative s-waves, although both pions are produced in a p- 
wave. Since no evidence of any significant p-wave contri- 
bution (except partially between the two emitted protons) 
is observed in the differential cross section, the dominating 
partial wave for two-pion production near threshold must 
be the 1 5'o partial wave in the entrance channel. However, 
for non-zero analyzing powers at least two interfering par- 
tial waves are necessary. Hence, from the Roper excitation 
and decay we do not expect sizable non-zero analyzing 
power values. Indeed, the calculations of Ref. [3], which 
include also the two-pion production via the AA system, 
predict vanishing analyzing powers. This result, however, 
should be taken with some caution, since the calculations 
have been carried out without taking into account initial 
state interactions, which should be a prime source of imag- 
inary parts of the amplitudes and analyzing powers are 
proportional to the imaginary part of interference terms. 

On the other hand analyzing powers originate from the 
interference of amplitudes. Hence small amplitudes in the 
production process, in particular p-wave amplitudes, can 
be easily sensed in polarized beam measurements 

2 Experiment 

The measurements have been carried out at the Jiilich 
Cooler Synchrotron COSY at beam energies of nominally 
747 and 793 MeV using the time-of-flight spectrometer 
TOF at an external beam line. The setup of the TOF de- 
tector system is displayed in Fig. 1. At the entrance of the 
detector system the beam - focussed to a diameter smaller 
than 2 mm - hits the LH2 target, which has a length of 
4 mm, a diameter of 6 mm and 0.9 /im thick hostaphan 
foils as entrance and exit windows |7j. At a distance of 
22 mm downstream of the target the two layers of the 
start detector (each consisting of 1 mm thick scintillators 
cut into 12 wedge-shaped sectors) were placed followed by 
a two-plane fiber hodoscope (96 fibers per plane, 2 mm 
thick fibers ) at a distance of 165 mm from target, see 
Fig. lb. Whereas the start detector mainly supplies the 
start times for the time-of-flight (TOF) measurements, the 
fiber hodoscope primarily provides a good angular resolu- 
tion for the detected particle tracks. In its central part the 
TOF-stop detector system consists of the so-called Quirl, 
a 3-layer scintillator system 1081 mm downstream of the 
target shown in Fig. lc and described in detail in Ref. |8J 
- and in its peripheral part of the so-called Ring, also a 3- 
layer scintillator system built in a design analogous to the 
Quirl, however, with inner and outer radii of 560 and 1540 
mm, respectively. Finally behind the Quirl a calorimeter 
(Fig. la,d) was installed for identification of charged par- 
ticles and of neutrons as well as for measuring the energy 




E cfllo [MeV] E cal0 [MeV] 

Fig. 2. AE — E plots for particles hitting the calorimeter. The 
AE information is taken from the TOF information obtained 
by Quirl and start detector and plotted as 1//3 2 , where P — v/c 
denotes the particle velocity relative to c. For the E informa- 
tion the energy deposited in the calorimeter has been taken. 
Left: without any trigger selections and conditions. Right: with 
trigger selection for three- and four-prong events as well as 
with missing mass conditions MM PP > 260 MeV and 100 MeV 
< MM ppvr < 180 MeV, see text. 



of charged particles. The calorimeter, details of which are 
given in Ref. [9] , consists of 84 hexagon-shaped scintillator 
blocks of length 450 mm, which suffices to stop deuterons, 
protons and pions of energies up to 400, 300 and 160 MeV, 
respectively. 

The particle identification is done with the standard 
Z\E-E technique using the energy signal from the calorime- 
ter together with the TOF-signal of the Quirl as shown in 
Fig. 2 and as described in more detail in Ref.|10|. where 
also sample spectra are shown. 

In order to distinguish 7r + and ir~ particles we have im- 
plemented the delayed-pulse technique into the calorime- 
ter setup. By using multihit TDCs the decay of a fi + 
resulting from the decay of a 7r + stopped in one of the 
calorimeter blocks can be registered by an afterpulse ar- 
riving in the fis range after the prompt pulse. In order 
to prevent false afterpulses to be registered, which stem 
from prompt pulses arising from follow-up events, we use 
the last layer of the Quirl as veto for afterpulses. This 
technique for ir + identification had previously been in- 
stalled at the forward detector of PROMICE/WASA and 
CELSIUS/WASA, respectively, where it has been shown 
to work reliably [lT|[l2 j. Fig. 3 shows a sample spectrum 
of registered delayed pulses at a proton beam of T p = 793 
MeV. The spectrum is very clean following the expected 
exponential decay pattern of ^ + in matter with a lifetime 
of 2.2 /is. The overall efficiency of this tt + identification has 
been determined to be 39 %. The efficiency losses mainly 
stem from up-stream moving positrons, which originate 
from /i + decay in the calorimeter and create a veto signal 
in the last layer of the Quirl. 

In the experiment the trigger for two-pion production 
was set to at least three hits in the Quirl-Ring system 
associated with two hits in the start detector. That way 
single-pion production events could be excluded, which at 
most can give two hits in the Quirl-Ring system. For cal- 
ibration and monitoring purposes, however, a two-prong 
trigger was also set up with a high prescale factor. From 
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straight-line fits to the hit detector elements the tracks of 
charged particles were reconstructed. They are accepted 
as good tracks, if they originate in the target and have 
a hit in each detector element the track passes through. 
In this way the angular resolution is better than 1° both 
in azimuthal and in polar angles. By construction of the 
calorimeter a particle may hit one or more calorimeter 
blocks. The number of blocks hit by a particular parti- 
cle is given by the track reconstruction. The total energy 
deposited by this particle in the calorimeter is then just 
the (calibrated) sum of energies deposited in all blocks 
belonging to the particular track. 

In order to have maximum angular coverage by the 
detector elements and to minimize the fraction of charged 
pions decaying in flight before reaching the stop detec- 
tors, the short version of the TOF spectrometer was used. 
In this way a total polar angle coverage of 3° < & < 
49° was achieved with the central calorimeter covering the 
region 3° < lab < 28°. 

By identifying and reconstructing three or four charged 
tracks of an event the reaction channel ppir + ir~ can be 
fully identified and reconstructed kinematically complete 
with 1-4 overconstraints. The maximum possible lab- 
oratory (lab) polar angle for protons emerging from two- 
pion production events is w 30°. Hence practically all pro- 
tons are within the angular acceptance of the calorimeter 
and can be identified as protons there. For charged pi- 
ons emerging from two-pion production events there is 
no kinematic limit at the incident beam energies of in- 
terest. Hence the angular coverage has not been complete 
for their detection with this setup. However, since detec- 
tion of two protons and one pion are enough to safely 
reconstruct the second pion of an event and since due to 
the identity of the initial collision partners the unpolar- 
ized angular distributions are symmetric about 90° in the 
overall center-of-mass system (cms), still most of the re- 
action phase space has been covered by this experiment. 
The total acceptance has been 15 % of the full phase space, 
which is a factor 5 -7 larger than in the PROMICE/ WAS A 
measurements, where the angular acceptance was < 21°. 

Though most of the events stemming from single-pion 
production could be removed by the three- and four-prong 
trigger conditions, there is still appreciable background 
left fromppir events, where the it either undergoes Dalitz 
decay or where the 7s emerging from ir° decay undergo 
conversion in matter. In both cases e + e~ pairs are pro- 
duced, which easily may be misidentified as charged pions. 
And since the ppir production cross section is three or- 
ders of magnitude larger than that for pp7r + 7r~ this back- 
ground is sizeable. In order to suppress this background we 
applied two constraints on missing mass spectra as shown 
in Fig. 4. For the missing mass MM vp7T of two identified 
protons and one identified pion, which has to be equiva- 
lent to the mass of the second pion in case of a true two- 
pion production event, we require 100 MeV < MM pp7T < 
180 MeV. We see from Fig. 4 that by this condition the 
MM pp spectrum gets already very clean. To be on the 
safe side we introduce as further constraint MM pp > 260 
MeV before applying the kinematic fit to events, who have 
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Fig. 3. Time spectrum of afterpulses detected by the delayed 
pulse technique installed at Quirl and Calorimeter. The log- 
arithmic slope given by the solid curve corresponds to r = 
2.18(6) /is, which agrees very well with the muon lifetime. In 
the fit a small background of 10 counts per bin (1 bin corre- 
sponds to 256 ns) is assumed. 
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Fig. 4. Left: missing mass spectrum MM pp7I of two identified 
protons and one identified pion with the condition MM VV > 
260 MeV. Right: missing mass spectrum MM PP of two iden- 
tified protons with the condition 100 MeV < MM ppn < 180 
MeV. 



passed both these conditions. For the kinematic fits we 
have one overconstraint in case of three-prong ppn events 
and 3 overconstraints in case of four-prong events, where 
all emerging particles have been recorded in the detector. 

The luminosity of the experiment was determined as 



dL/dt = 5 * 10 



2* 



cm 



_1 from the analysis of pp elastic 



scattering, which was measured simultaneously utilizing a 
prescaled two-prong trigger for hits in the Ring detector. 



Due to their opening angle of S p 



80° between both 



tracks, such two-track events have both hits in the Ring. 
They are easily identified by using in addition the copla- 
narity constraint 170° < A<P < 180°, i.e. using the same 
procedure as described in more detail in Ref. [10]. Ad- 
justment of the pp-elastic data in absolute height to the 
SAID values gives the required luminosity. All data have 
been efficiency corrected by MC simulations of the de- 
tector setup by using the CERN GEANT3 Q3] detector 
simulation package, which accounts both for electromag- 
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netic and hadronic interactions of the ejectiles with the 
detector materials. 

The measurements have been conducted with a vector- 
polarized proton beam providing a steady change of "spin- 
up" and "spin-down" runs in subsequent spills of about 
10 s duration. "Spin-up" is defined of having the spin 
of the protons perpendicular to the plane defined by the 
COSY accelerator ring and pointing to the sky (y-axis). 
According to the Madison convention [14] we then use 
a coordinate system, where the z-axis is along the beam 
direction and the x-axis points to the horizontal left cor- 
responding to the azimuthal angle of <P = 0° being at 9 
o'clock when looking downstream the beam direction. 

For the determination of the beam polarisation at the 
target we use again elastic scattering, determine its coun- 
trate asymmetries and compare them to the known ana- 
lyzing powers of elastic scattering as compiled in the SAID 
database [TS]. The countrate asymmetry e(0,<P) at a spe- 
cific solid angle (<9, <P) for a certain particle of an event is 
given by the difference of countrates at " spin-up "7V|(0, <&) 
and "spin-down" Ni(0,<I>) normalized to their sum: 



e(0,<P) 



_ jv T (e,'g)-jv^(e,'g) 



= Py*Ay(0,-P), 



(1) 



where P y is the beam polarisation and 

A y (0,$)=A y (0)*cos($) (2) 

is the analyzing power of the reaction of interest with the 
particle of interest detected at the solid angle (0, <P). 

In order to determine and also to monitor online the 
beam polarisation we have selected pp elastic scattering 
events with 32° < iah < 33°, where the SAID database 
[15J gives A y = 0.30. A sample of the measured countrate 
asymmetry e(<9, <P) is given in Fig. 5 exhibiting the ex- 
pected cosine dependence on the azimuthal angle. A fit to 
the data with eqs.(l) and (2) yields the beam polarisation 
P y . During the beamtime the beam polarisation was very 
stable within a few percents and with very high values of 
typically P y = 0.87. 

Due to its cylinder symmetry the TOF detector is ide- 
ally suited for polarisation measurements. Thus, the an- 
alyzing power A y may be determined experimentally by 
use of the super ratio R{0, ( P) of the four countrates for 
the combinations of "spin-up", spin-down" with "left" or 
"right" : 



R(0,<P) 



N r (0,<P)*N l (0 : <j>+18O° ) 
N l (0,<P)*N r (6>,<J>+180° ) 



and 



A (0 $\ - J- * R ( e ^')- 1 



(3) 



(4) 



In this superratio systematic uncertainties, detector ef- 
ficiencies etc. cancel to a very large extent making thus 
the deduced analyzing power essentially free of systematic 
asymmetries, which often are not under control in "single 
arm" measurements with a not cylinder-symmetric detec- 
tor. By using the four countrates we in addition may define 
an "anti-analyzing power" (instrumental asymmetry) A y 
with the alternative super ratio R defined as 
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Fig. 5. Countrate asymmetry e(<9i a t>,$) in dependence of the 
azimuthal angle for elastic scattering events measured in the 
polar angle interval 32° < Qi a b < 33°. The solid curve is a fit 
to the data with a cos(&) dependence. 
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Fig. 6. Analyzing power of elastic scattering events measured 
in the Ring detector in comparison with the SAID database 
[15] (solid line). 



which leads to A y = for all cases with the exception of 
systematic changes in beam direction when flipping the 
spin of the projectiles. Hence A y serves as a control ob- 
servable to uncover systematic asymmetries of such kind. 

Due to their trivial cos^ dependence the analyzing 
powers are usually quoted as 

Ay(0) = Ay(0,^ = 0). 

In order to take into account the fully available statis- 
tics, we have sorted the data into and bins, calcu- 
lated A y (0,<l>) for each of these bins and fitted for each 
particular bin the <P dependence of A y (0,<P) by a co- 
sine dependence. As a result we obtain for each bin 
an A y (0) value, which contains the full statistics of the 
measurements over the full <P range. 

As an example we show in Figs. 6 and 8 analyzing 
power and anti-analyzing power (instrumnetal asymme- 
try) distributions obtained for pp elastic scattering events 
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Fig. 7. Analyzing power data for the pp — * dir + reaction as 
measured in the TOF detector (solid circles) in comparison 
with the SAID database [E] (solid line) and LAMPF data [E] 
(asterisk symbols). 
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Fig. 8. "Anti-analyzing power" (instrumental asymmetry) 
as denned in eq.(5) for elastic scattering (left) and two-pion 
production events in dependence of the two-pion polar angle 

or (right). 



3 Results 

3.1 Cross Sections - Unpolarized 

Due to the identity of the collision partners in the entrance 
channel the angular distributions in the overall center-of- 
mass system have to be symmetric about 90 °, i.e. the full 
information about the reaction channels is contained al- 
ready in the interval 0° < cm < 90° . Deviations from this 
symmetry in the data indicate systematic uncertainties in 
the measurements. Hence we plot - where appropriate - 
the full angular range, in order to show the absence of 
such type of systematic errors in our measurement. 



Tab. 1: Total cross sections Otot at T p « 750 and 800 
MeV for the reaction pp — » ppn + ir~ evaluated in this work 
and compared to previous measurements. 









T p [MeV] 


this work 


previous 


750 


1.6 (2) 


1.6 (3) Q ) 






1.3 (3) 6 > 


800 


4.1 (4) 


3 (l) c ) 


Q ) Ref. [H] 


b) Ref. [§] c ) Ref. 


ED 



The evaluated total cross sections for the ppir + ir~ chan- 
nel are given in Tab. 1 and shown in Fig. 9 together with 
previous results. The uncertainties assigned are based on 
systematics for acceptance and efficiency corrections as 
obtained by variation of MC simulations for the detector 
response, where we have varied the MC input assuming 
either pure phase space or some reasonable models for the 
reaction under consideration. For 750 (800) MeV we have 
collected about 1600 (6000) good events, i.e., statistical 
uncertainties are of minor importance for the uncertain- 
ties in the absolute total cross section. 



recorded in the Ring detector. The analyzing power data 
agree well with SAID in the measured range. As a further 
check of the deduced beam polarisation we have evaluated 
also the pp — » dn + reaction utilising the same 2-prong trig- 
ger as used for the selection of elastic scattering events. 
The results for the selected d7r + events is shown in Fig. 7 
in comparison with the SAID data base [15] and LAMPF 
data [16]. Again there is good agreement. 

The anti-analyzing power (intsrumental asymmetry) 
data are all consistent with zero giving thus no indication 
of systematic asymmetries in the measurements. As a typ- 
ical example of an anti-analyzing power distribution ob- 
tained from two-pion production events we show in Fig. 8, 
right, the distribution according to the angle <9£" 1 of any of 
the pions in the overall cms. Again the data are consistent 
with zero, whereas A y {0™ 1 ) / 0, see Fig. 15. 



For T p = 750 MeV total and differential cross sections 
for this reaction channel are very well known from pre- 
vious experiments at CELSIUS [4,18]. Hence we use the 
analysis of our data at this energy primarily as a check of 
the reliability of our measurement and data analysis. In 
fact, the differential distributions obtained at CELSIUS 
and COSY-TOF are compatible within statistical uncer- 
tainties [9,19j. As an example we compare in Fig. 10 both 
data sets for the differential distributions of the invariant 
mass Man and the opening angle 6^ between the two pi- 
ons in the overall cms. We also note that the values for 
a tot agree very well within uncertainties. 

In contrast to T p = 750 MeV there exists merely a 
single-arm measurement for T p = 800 MeV, namely the 
one performed at LAMPF [if] with an estimated total 
cross section of 3 /ib. 

Differential cross sections for invariant masses and an- 
gular distributions obtained from our measurements are 
shown in Figs. 10 - 13. 
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Fig. 9. Energy dependence of the total cross section. Shown 
are the data of this work (solid squares) together with those 
from CELSIUS [4ll5ll6] (solid circles), Gatchina - old [23] (open 
circles) and new [24] (solid triangles), LAMPF [TTpCi] (open 
diamond and star, respectively), KEK [22] (open squares) and 
Berkeley [21] (open triangles). Solid and dashed curves corre- 
spond to theoretical calculations of Ref.[3] with and without pp 
FSI. The hatched area represents the phasespace dependence 
adjusted arbitrarily to the cross section at T p = 750 MeV. 
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Fig. 10. Differential cross sections for the distributions of tttt 
invariant mass Mm (left) and mv opening angle in the overall 
cms (right) at T p = 750 MeV. The data of this work (filled 
symbols) are compared to those obtained at CELSIUS (open 
symbols) as well as to phase space distributions (shaded areas) 
and calculations according to eq.(6) (solid lines). CELSIUS 
data and theoretical curves are renormalized in area (utot) to 
the data of this work. 



3.2 Analyzing Powers 
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Fig. 11. Results of this work for differential distributions of the 
invariant masses M pp ,M nn ,M p7T and M p7r7r at T p = 800 MeV. 
They are compared to phasespace distribuions (shaded areas) 
as well as to calculations according to eq.(6) (solid lines). 
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Fig. 12. The same as Fig. 11, but for of the opening angle 
between the two pions Stttt, the polar angles of protons O p m , 
pions (9£ m and of the two-pion system (9£™ - all in the overall 
cms. 



These measurements are the first ones with polarized beam, 
which have been conducted for the two-pion production 
channel in NN collisions. The angular distributions for 
the analyzing power have been obtained as described in 
detail in section 2. 
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Fig. 13. Differential cross sections for the distributions of tt7v 
invariant mass (left) and irn opening angle in the over- 

all cms 8^^ (right) at T p = 800 MeV. The data of this work 
are compared to phase space distributions (shaded areas) and 
calculations according to Ref.[3j. The dashed lines show the 
original calculations of Ref. [3j renormalized to the data in 
area, the solid lines the result, when we reduce the amplitude 
for the branch N* — ► Air by a factor of two relative to the 
branch N* — + Na. 
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For T p = 750 MeV, where the statistics was quite mod- 
erate, we show in Fig. 13 our results for the A y distribu- 
tions on the angle of any proton 0£ m , of any pion <9£ m , 
of the 7T7T system and of any pir system (9^™ - all in 
the overall cms. 

For T p = 800 MeV, where we have much better statis- 
tics, we show in addition to these distributions (Fig. 15) 
also the ones, where we differentiated between n + and n~ 
particles with help of the delayed pulse technique (Fig. 
16). Finally we display in Fig. 17 the analyzing powers for 
the angles 0£ w of any pion in the irir subsystem and for 
the angles 0P p of any proton in the pp subsystem, both 
taken in the Jackson frame, i.e., having the beam axis as 
z-axis. 



4 Discussion of Results 
4.1 Cross Sections - Unpolarized 

The total cross sections obtained in this work are shown 
in Fig. 9 together with data from previous work [17,20, 
[2T1l22jl23p 24 [ l4) and theoretical calculations of Ref [3]. At 
T p = 750 MeV our value agrees well with those obtained 
at CELSIUS OH]. At T p = 800 MeV our value is consis- 
tent with the estimate of Ref. [17] and follows also the 
trend given by the very recent bubble chamber results 
from Gatchina [24], which correct the old Gatchina results 
|23] for the total cross sections. 

The theoretical calculations of Ref. [3] are compatible 
with the data from LAMPF [T7], CELSIUS gE], Gatchina 
[24] and COSY-TOF for T p < 800 MeV only, if set I (see 
Ref. [2]) for the Roper resonance parameters is used - and 
if the pp final state interaction (FSI) is neglected (lower 
dashed curve in Fig. 9) . Taking into account the FSI (solid 
curves) the theoretical calculations overestimate the cross 
sections in the near-to-threshold region by as much as an 
order of magnitude. 



Fig. 14. Analyzing power distributions for the polar angles 
&; m , 6>r, &l™ and 6>™ at T p = 750 MeV. 
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Fig. 15. Analyzing power distributions for the polar angles 
®p m , 6>£ m , ©™ and & c ™ at T p = 800 MeV. The solid curves 
show a fit according to eq.(7). 



For the discussion in particular of the data taken at 
T p = 800 MeV, we compare our data in Figs. 10 - 13 with 
phase phase distributions denoted by the hatched areas in 
the figures, with theoretical calculations of Ref. 0] (Fig. 
13 only) as well as with simplified calculations of Roper 
excitation and decay according to the ansatz for the Roper 
amplitude 



A - 1 + cki * k 2 (3D A ++ + D A o) 



(6) 
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Fig. 16. Analyzing power distributions for the polar angles 
6>X 6>l m , C ™ + and 6>;™_ at T p = 800 MeV, i.e. separated 
for positively and negatively charged pions. The curves show a 
fit according to eq.(7). 
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The coefficient c contains the branching of the decays 
N* -> Ncr and N* -» An. Adjustment of c to the CEL- 
SIUS data at T p = 750 and 775 MeV [5] resulted in a 
branching ratio of R = 1/^*— >ATt^>Nmr/ rN*—>Na between 
the two-pion decays of the Roper resonance of R = 3(1), 
if the pole of the Roper resonance is assumed to be at 
its nominal value of 1440 MeV - in good agreement with 
the PDG value of 4(2) [25], but more precise. However, if 
the pole is taken to be at 1371 MeV, which meanwhile is 
considered to be more correct [25,26,27,28,29], then this 
ratio changes to R = 1.0 (1). Calculations with this value 
of c are shown in Figs. 10 - 12 by the solid lines. They 
give a qualitative description of the main features in the 
differential data. At a closer look, however, they provide 
a shift in the M V7r spectrum and a slope in the Sn-n spec- 
trum, which is somewhat too large. This suggests that the 
A^* — > An decay is even weaker than the N* — * Na decay 
- as was also obtained in a very recent analysis of photon- 
and pion- induced reactions [29]. The result obtained here 
is also in agreement with a recent analysis of pp — > pp7r°7r° 
data taken at CELSIUS- WAS A [33]. 

Fig. 13 shows the full calculations of Ref.[3], which also 
include reaction processes other than just the Roper ex- 
citation. The dotted line corresponds to calculations with 
the same branching for the N* — > Na and N* — > Air de- 
cay processes as used with the simplified ansatz (6). We 
see by comparison of Figs. 11 and 12 with Fig. 13 that 
full and simplified calculations give nearly the same re- 
sult, which means that the terms neglected in the ansatz 
(6) are of minor importance. In fact, if we now reduce 
the N* — > Air amplitude by a factor of two in the calcul- 
tions of Ref.[3j - as done in Ref. [33] - then we get full 
agreement with the data. 



4.2 Analyzing Powers 



cos 6" 



Fig. 17. Analyzing power distributions for the polar angles 
0p P , OZ n in the pp and irir subsystem (Jackson frame) at T p = 
800 MeV. The curves show a fit according to eqs.(7) and (8), 
respectively. 



which multiplies the expressions for a exchange, FSI 
and N* propagator (see Refs. [llll]) and where Da denotes 
the A propagator. The first term in eq. (6) stands for the 
Roper decay into the Ncr channel. The second term gives 
the Roper decay via the A resonance, i.e. N* — » An — > 
Nnn, where the scalar product of the pion momenta ki 
and k 2 stands for the double p-wave emission of the pions 
in this process. Whereas the first term provides a phase- 
space like behavior of the observables, the second term 
is proportional to the cosine of the 7T7t opening angle in 
the cms and hence provides both a change in the slope 
of the spectrum and a shift of the M TI distribution 
towards higher masses. From this apparent signature in 
these observables the parameter c can be fixed easily by 
corresponding data. 



The analyzing powers taken at T p = 750 MeV shown in 
Fig. 14 are essentially compatible with zero within their 
large statistical uncertainties. Because of that we do not 
discuss them here in more detail and concentrate on the 
data at T p = 800 MeV, which have much better statistics. 
We only note in passing, that the trend in the analyzing 
powers at T p = 750 MeV is just opposite in sign to that in 
the data at T p = 800 MeV. The reason for it is not known 
to us. 

At T p = 800 MeV the data displayed in Figs. 15 - 
17 exhibit significant analyzing powers with values up to 
A y = 0.3. Whereas the analyzing powers in the pp and 7T7t 
subsystems are very small and compatible with zero, they 
are quite substantial in the overall cms for (9£ m and <9£™ . 
The analyzing power of the protons in the overall cms, i.e. 
A y (0 p m ), is again small. 

The analyzing powers for (9£ m and 6>£™ are of course 
not decoupled. Since for the latter the analyzing power 
is largest, the source for polarization producing effects is 
expected to be connected with this observable. 

From general considerations [30,31] of meson produc- 
tion in NN collisons the analyzing powers in the overall 
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cms has the following dependences (see appendix for 
more details): 

A y = a sin(0) + b sin(20), (7) 

whereas in the pp subsystem we have only a 20 depen- 
dence: 

A y = V sin{20). (8) 

A more detailed description within the partial wave con- 
cept is given in the appendix. The solid lines in Figs. 15 - 
17 show fits to the data with eqs. (7) and (8), respectively. 
As one can see from Figs. 15 - 17, the coefficients a and b 
are only significantly different from zero for 0p m , 0% m and 
(9£™ distributions. In addition the latter ones are essen- 
tially described by a sin(0) dependence with the sin(20) 
dependence being of minor importance or negligible. 

Finally we look into possible differences between 7r + 
and 7T~ angular distributions by exploiting the delayed 
pulse technique to identify the 7r + particle in an event. 
That way we obtain the data distinguished between 7r + 
and 7r~ observables in Fig. 16. Though the statistics gets 
worse due to the finite tt + identification efficiency, we find 
the analyzing powers for the ir + observables in 6>^™ and 
0p™+ to be systematically larger than for those for the 
corresponding ir~ observables. Such a difference indicates 
a significant contribution of isovector irir pairs. 

From this and the requested sine dependence of the 
analyzing power we learn that the driving term for finite 
analyzing powers must be connected primarily with the 
7T + momentum and be of the form 

a- (q X (9) 

where q and k ff + are momentum transfer between the nu- 
cleons and the ir + momentum, respectively, and a is the 
Pauli spin operator. Since the momentum transfer close to 
threshold is approximately given by the initial momentum 
in the overall cms, we arrive with such an operator at the 
observed sin((9£™) dependence. 

Since to lowest order spinflip terms like (9) are con- 
nected with p-waves, whereas nonspinflip terms are con- 
nected with s-waves, the analyzing power arises to lowest 
order from an interference of S- and P-waves in the en- 
trance channel. The dominating reaction process at T p « 
800 MeV is via the Roper excitation and hence fed by the 
1 Sq partial wave in the entrance channel. Processes involv- 
ing the operator of eq. (9) evolve via the 3 Pi partial wave. 
From inspection of the formulas given in Ref. [3] we see 
that neither Roper nor AA excitation and decays provide 
such an operator. The only process, which provides such 
an operator in the theoretical investigations of Ref. [3], is 
single A excitation and decay in combination with pion 
rescattering on the other vertex as given by the graphs 
(10) and (11), Fig.l of Ref. [3], or graph (9), which gives a 
successive excitation of a A on a single nucleon . However, 
in these calculations such graphs give only tiny contribu- 
tions to the cross section. On the other hand the analyzing 
power is an ideal tool to reveal such small amplitudes by 
their interference with large amplitudes. Hence it appears 
very interesting to check by realistic calculations, whether 
such graphs could be the origin of the finite analyzing 



powers we observe. As pointed out above, the calculations 
of Ref. [3] are done in the plane wave limit, i.e., are not 
valid for polarization observables. Hence an improved cal- 
culation, which includes initial state interactions, would 
be of great help to understand the polarization data on a 
quantitative basis. 

Another interesting aspect could be the consideration 
of a possible influence of the Z\(1600) resonance. Though 
its resonance energy is as high as that of the AA excita- 
tion, its unusually large width of « 350 MeV [25] is in favor 
of taking part in the two-pion production already close to 
threshold. Indeed, as has been noted already in Ref. |32j, 
the isotensor two-pion production amplitude is observed 
to be unusually large. And the Z\(1600) resonance is the 
only resonance candidate to contribute to the isotensor 
amplitude near threshold in addition to the AA process. 

The analyzing power arises from the interference of 
spinflip and non-spinflip terms in the reaction amplitude, 
whereas the unpolarized cross section is just the sum of 
their moduli squared. If we denote the non-spinflip term 
by g and the strength of the spinflip term which multiplies 
the operator (9) by h, then we have 

a A y = -2Im{g*h) fc 7r + sin(&%?) (10) 

and 

a = | 5 | 2 + H 2 fc2 + sin 2 (9^) < \g\\l + 4/4), (11) 

i.e., a 30% effect in the analyzing power produces only 
a 2% effect of non-isotropy in the differential cross section. 
Such an effect is not observable in the cross section with 
the statistics accumulated in this measurement. 



5 Summary 

The experimental results presented in this work consti- 
tute the first exclusive measurements of 7r + 7r~ production 
in NN collisions at T p = 800 MeV and provide the first 
polarization data at all for this production process. The 
unpolarized differential cross sections support the findings 
of previous analyses that Roper excitation and decay is the 
dominant two-pion production process. The data favor the 
direct decay of the Roper resonance into the Na channel 
over its decay via the A resonance. This in turn supports 
the breathing mode nature of the Roper resonance. 

The observation of sizable analyzing powers in com- 
bination with the fact that neither Roper nor AA pro- 
cesses contribute to substantial analyzing powers, point 
to a small amplitude of different nature, which interferes 
with the strong Roper amplitude. Candidates for such a 
small amplitude could be A excitation processes, where 
on the one nucleon vertex a Z\(1232) orZ\(1600) is excited 
followed by its transition A — > An — > Nnn, which would 
be the first time that such a process has been sensed in 
NN collisions. 

This work has been supported by BMBF, DFG (Europ. Gra- 
duiertenkolleg 683) and COSY-FFE. We acknowledge valuable 
discussions with Eulogio Oset, Murat Kaskulov and Alexander 
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Fig. 18. Definition of the coordinate system. Dashed (solid) 
lines denote pions (nucleons). 



Sibirtsev. We are particularity indebted to Luis Alvarez-Ruso, 
who made his computer code available to us. 



6 Appendix 

In this appendix the general structure of the amplitudes 
is derived. Since none of the observables shows a stric- 
ing energy dependence and the energy is still not too far 
above threshold it appears justified to work with the low- 
est partial waves only. Hence we will only derive the an- 
gular structures of differential cross section and analyzing 
power up to order Q 2 , were Q is any of the momenta in 
the final state as defined in Fig. [TH 

Before we can identify the partial waves that are rele- 
vant given the assumption formulated above, we need to 
briefly comment on the relevant selection rules. The Pauli 
principle enforces that a proton pair of spin 0(1) has even 
(odd) parity. As a consequence, for unpolarized observ- 
ables pp S- waves and pp P-waves never interfere. Since 
in the experiment presented the polarization of the outgo- 
ing protons is not measured, no interferences of pp P- and 
S-waves in the final state are possible. We therefore only 
need to consider pp P-waves in the final state in combina- 
tion with s-waves for the other two subsystems (7T7t and 
7T7T relative to pp) . 

There are 4 partial waves possible with an NN p-wave 
in the final state ( 3 Pj — > 3 P/s[s] with j = 0,1,2 and 
3 P2 — > 3 P2s[s], where we used the notation 2S+1 Lj — » 
2S +1 L'j,T[l], with S, L, j (S", L' , j) for spin, angular mo- 
mentum, and total angular momentum of the pp pair in 
the initial (final) state, I for the angular momentum of the 
pion pair, and / for the angular momentum of the pion pair 
relative to the pp pair. It is straight forward to see that 
only 3 independent operator structures are possible and 
we define 

M p = d 1 (S-p)(S'-p')+d 2 (S-p')(S'-p)+d 3 (S-S'){p-p') . 

The partial wave amplitudes are linear combinations of 
the di. 

Terms of at most quadratic dependence on the final 
momenta in combination with pp 5-waves in the final 
state can emerge from overall 5"-waves interfering either 
with themselves, p-waves in the subsystems or d-waves in 



the subsystems. We therefore need to consider in total 6 
partial waves for this part that we parametrize as 



M S = a + h(S X p) ■ q' + b 2 ( (q' ■ pf - \q' 2 



+ c 1 (Sxp)-k' + c 2 (q'-k') + 
+ c 3 Up-q')(p-k')-~( q '-k')) , 

where the a term corresponds to 1 Sq — > 1 5qs[s], the b\ 
term to 3 Pi — > 1 Sop[s], the b 2 term to 1 D 2 — > 1 S , d[s], the 
ci term to 3 Pi — > 1 5os[p], the c 2 term to 1 Sq — > 1 Sop[p], 
and the C3 term to 1 D 2 — > 1 Sop\p]. The full matrix element 
is given by 

M=M S + M P . 

Using the formulas of chapter 4.3 of Ref. (31] it is 
straightforward to calculate both unpolarized cross sec- 
tions as well as analyzing powers from the expressions 
given above. We find (dropping terms more than quadratic 
with respect to the final momenta) 

4a = \a\ 2 + Ib^q' 2 sm(9 ql ) 2 + \ Cl \ 2 k' 2 sm(9 k ,) 2 + 

+ 2Re{a*b 2 q' 2 (cos(6 q ,) 2 - 1 
2 

+ -a*c 3 q'k' cos(6 q i)cos(6 k >) + 

a c 2 - -a c 3 + b x c\ I k q 
sm(0fc/)sin(0g')cos((^/ -</>»,/)} 

+ P ' 2 (M 1 | 2 + |d 2 | 2 ) + 

+ p' 2 cos(^) 2 (3Ma| 2 + 2ReKd 2 + d\d z + d* 2 d 3 }) 
4<j Ay = - Im {{d\d 2 + d* ± d 3 - d 2 d 3 )p' 2 sin(20 p /) cos(0 p /) 

+2b\aq sm{6 q i) cos(0 g <) + 2c*afc' sin(0fc') cos(0fe')} 

Note that, as a consequence of the fact that two ampli- 
tudes with a spin singlet entrance channel can not interfere 
in the analyzing power, the expression for the polarization 
observable has a much simpler structure. 
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